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Abstract 

We study the total gravitating mass distribution in the central region of 23 clus- 
ters of galaxies with Chandra. Using a new deprojection technique, we measure the 
temperature and gas density in the very central region of the clusters as a function 
of radius without assuming any particular models. Under the assumptions of hydro- 
static equilibrium and spherical symmetry, we obtain the deprojected mass profiles 
of these clusters. 

The mass profiles are nicely scalable with a characteristic radius (r 2 oo) and mass 
(M200) on the large scale of r > 0.1r 2 oo- In contrast, the central (r < 0.1r 2 oo) mass 
profiles have a large scatter even after the scaling. The inner slope a of the total mass 
density profile (p(r) oc r~ a ) is derived from the slope of the integrated mass profile. 
The values of the inner slope a at the radius of 0.02r2oo («o) span a wide range from 
to 1.2. For 6 out of 20 clusters, «o is lower than unity at a 90 % confidence level. 
CDM simulations predict that the inner slope a is in the range 1 < a < 2, which is 
inconsistent with our results. We also found that the gas fraction near the center of 
a cluster has a negative correlation with a - Our result suggests that the gas-rich 
clusters in the central region tend to have a flat core. 

Key words: galaxies: clusters: general - cosmology: dark matter - X -rays: galax- 
ies: clusters 
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1. Introduction 



The Cold Dark Matter (CDM) model has become the standard paradigm for explaining 
observations of the large-scale structure of the universe. In the CDM model, dark matter con- 
sists of non-baryonic, collisionless, cold particle. The properties of dark matter density profiles 
in the CDM model have been investigated extensively through numerous N-body simulations. 
Navarro, Frenk, & White (1997) (hereafter NFW) claimed that the dark matter density profiles 
in the CDM model are reasonably approximated by a universal form with singular behavior in 
its central region. Several N-body simulations predict that the density of dark matter increases 
as a power law p(r) oc r~ a , with a in the range of 1 to 2, in the central region (e.g., a = 1 
by NFW; a = 1.5 by Moore et al. (1998)). Measurements of the inner slope a of dark matter 
density profiles offer a powerful test of the CDM model. 

The observational efforts in this respect have been in the form of dynamical studies of 
low surface brightness and dwarf galaxies. The observations obtained in those studies suggest 
the presence of a relatively flat core: < a < 1 (e.g., Firmani et al. (2001)). Gravitational 
lensing has made some observational constraints available at the scale of galaxy clusters. For 
instance, Sand, Treu, & Ellis (2002) showed that steep inner slopes (a > 1) are ruled out at 
better than 99 %, for the lensing cluster MS2137-23. Although gravitational lensing studies 
provide a unique and important probe of dark matter profiles, they generally can be applied 
only to a limited sample of clusters that satisfy a specific lensing condition. X-ray observations 
of the density and temperature of a hot intracluster medium (ICM), on the other hand, probe 
the mass of a cluster of galaxies under the assumption of hydrostatic equilibrium. This could be 
a powerful tool to investigate dark matter profiles in the central regions of clusters. However, for 
previous X-ray satellites, such as ROSAT and ASCA, the detailed study of ICM temperature 
and density profiles at small scales has been difficult because of limitations on the performance 
of imaging or spectroscopic instruments. The high spatial resolution imaging spectroscopy of 
Chandra enables the measurement of mass profiles in the very central regions of clusters of 
galaxies. Several groups have obtained X-ray constraints on the dark matter profiles of some 
clusters. These results are apparently consistent with the CDM model (e.g., David et al. (2001), 
and Arabadjis, Bautz, & Garmire (2002)). However, Ettori, Fabian, Allen, & Johnstone (2002) 
showed that the mass profile of A1795 flattens within 100 kpc. Systematic studies are thus 
required for a large sample of clusters. In this paper, we systematically study the mass profiles 
of 23 clusters of galaxies. 

We assume Q m = 1, Q\ = 0, and H = 50 km s" 1 Mpc" 1 throughout this paper. Unless 
otherwise noted, all errors are la (68.3 %) confidence intervals. 
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2. Sample and Analysis 



We selected our sample from Chandra archival data of galaxy clusters. To obtain spa- 
tially resolved spectra, we restricted the observations to those in which ACIS were employed 
without gratings. By the end of September 2002, the archive contained observation data for 
about 150 clusters (~ 200 pointings) that met this criterion. 

We applied the following criteria for further selection of the data in order to meet our 
main concern, investigation of the central mass profiles of galaxy clusters. First, clusters must 
be bright enough to provide spatially resolved spectra with good statistics. For this we referred 
to the catalog of Reiprich & Bohringer (2002), which consists of 106 bright clusters compiled 
from several catalogs based on the ROSAT All-Sky X-ray Survey (Voges et al. (1999)). The 
minimum X-ray flux among the Reiprich samples is 0.234 x 10 -11 ergs s _1 cm -2 (0.1-2.4 keV). 
This is bright enough for our analysis under typical observational conditions. Among the 106 
clusters in Reiprich sample, 43 clusters are included in the Chandra data archive. Secondly, 
clusters should be spherically symmetric as our deprojection analysis assumes the spherical 
symmetry. We thus excluded merging clusters like A754 (Henriksen & Markevitch (1996)). 
Although this second criterion is somewhat ambiguous, we will examine how this spherical 
symmetry assumption affects the final result in Section 5. The last criterion is that the X- 
ray emission from the outer region of a cluster must be covered by the detectors used in the 
observations. This is because the deprojection analysis depends on accurate measurement of 
temperatures and densities of the outer regions of clusters. Data for 20 clusters met all of these 
criteria. We also employed three distant clusters, A1835, A963, and ZW3146, that are bright 
and spherically symmetric but that are not included in Reiprich & Bohringer (2002). 

The observation log and the properties of each cluster are summarized in Table 1 and 
Table 2, respectively. The redshifts of the 23 sample clusters range from 0.0110 to 0.2906, with 
a median of 0.0852. 

Data reduction and analysis were performed with the Chandra Interactive Analysis of 
Observations package, CIAO-2.2, with calibration database CALDB-2.12, as provided by the 
Chandra X-ray Center (CXC). We started the reduction from the standard level 2 event files 
archived at CXC, which are the products of the pipeline processing. We adopted the standard 
reduction scheme by following the CI AO threads 1 . 

To remove the flare events, we performed lightcurve screening using the CIAO task 
lc-dean. We made a background lightcurve, a time history of the event rate taken from a source- 
free region on the detector, with a time bin size of 259.28 s using the CIAO task lightcurve. In 
order to exclude flare events, we discarded the data taken at the time the count rate deviates 
from the mean by ±3<r, where the mean value is defined during the quiescent period. Point 
source detection was performed with the CIAO wavelet source detection routine wavdetect with 

1 http: / / asc.harvard.edu/ciao/threads 
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a significance parameter 10 . We made a 0.3-10 keV image binned by using a bin size of about 
2x2" (4x4 pixels). The area around the detected point sources was excluded in the following 
analysis. In order to estimate the background level to be subtracted from the X-ray spectra and 
images, we applied the blank-sky data compiled by Markevitch (2001) 2 as background data. 
These background data are event files made with the same lightcurve screening process of the 
cluster data. 

The spectra were extracted in the concentric annuli centered on the X-ray peak with 
different widths to ensure similar statistics in the background-subtracted spectra. The X- 
ray peak was determined with the X-ray images from which point sources were removed. We 
examine what the appropriate setup is for the width of the annuli, or equivalently, the statistics 
of each spectrum, using a simulation. From this simulation, we found that the lower photon 
counts or the more annulus result in the larger systematic errors in the temperature profile as 
shown in Arabadjis, Bautz, & Garmire (2002). In order to suppress the systematic errors in the 
temperature profile, we restricted the photon count per each annulus and the number of annuli, 
as follows: (1) The photon count per each annulus must be at least 1 x 10 4 , and (2) The number 
of annuli N must be 5 < N < 10. The radius of the outermost annulus was determined to cover 
the 4-<t background level of the ROSAT PSPC image. The average of the outermost radius 
is about 720 kpc. The background spectra were extracted from the background data with the 
same regions on the detector. Redistribution Matrix Files (RMF) and Auxiliary Response Files 
(ARF) were made using the CIAO tasks of mkrmf and mkwarf. These tasks make a weighted 
RMF and a weighted ARF for the spectral analysis based on a 32 x 32 pixel grid of calibration 
files. This is because the RMF and ARF vary with detector location. To compensate for the 
degradation in low-energy efficiency, we used the tool corrarf provided by CXC. The corrarf 
corrects the ARF according to the observation date. 

2.1. Deprojection Analysis 

To determine the deprojected temperature and gas density profiles, we applied a new 
deprojection technique developed by Arabadjis, Bautz, & Garmire (2002). We here briefly 
summarize this technique (see also Arabadjis, Bautz, & Garmire (2002) and Katayama (2003)). 

A Schematic view of the deprojection analysis is shown in Figure 1 (Left). In this 
example, we extract spectra from N concentric annular regions. The projected luminosity Sj 
in a given energy band on the jth annulus is expressed by the integration of emissivities along 
the line of sight. The relationship between Sj and the volume emissivity of the ith spherical 
shell is expressed as 

^ = XX-e, (1) 

i=j 

where V^- is the volume of the ith spherical shell intersected by a cylindrical shell whose radius 
2 http://cx c. harvard.edu/contrib/maxim/bg/index. html 
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equals the jth projected annulus. Note that we have to make sure, or make the assumption that, 
X-ray emission is negligible outside of the outermost annulus. Vij is geometrically calculated as 

Va = \v[{r 2 i+1 - 6?) 3/2 - (rl, - b] +l fl> - (rf - b))^ + (r 2 - b) +l )^] (i > j) 

= (i<j), (2) 

where and r i+1 are the inner and outer radii of the ith spherical shell, and bj and bj + i are the 
inner and outer radii of the jth annulus, which equal rj and 7j-+i, respectively. Since Equation 
1 can also be written for all annuli, these are written as 

S = V-e. (3) 

Since V is a triangle matrix, we can obtain the e by solving the inverse matrix V -1 . 

In most previous deprojection analyses, only X-ray spatial information was utilized and 
an additional assumption on the temperature profile kT(r) or the potential profile <f>(r) was nec- 
essary. Even when both types of information are available, some authors assume the potential 
profile 0(r), or equivalently, the gravitational mass density profile p(r), beforehand. However, 
the method by Arabadjis, Bautz, & Garmire (2002) does not make such assumptions. We first 
make a trial model for volume emissivity at each spherical radius, which is a function of gas 
temperature fcTj, gas density n g ^ and gas abundance Z» when we employ an X-ray emissivity 
model of thin thermal plasma. We adopted the MEKAL (Mewe, Gronenschild, & van den 
Oord 1985; Mewe, Lemen, & van den Oord 1986; Kaastra & Mewe 1993; Liedahl, Osterheld, 
& Goldstein 1995) model in the XSPEC data analysis package for our X-ray emissivity model, 
in which normalization Ki is used instead of gas density n g i . Therefore, the number of free 
parameters to be determined is 3xN except for an additional free parameter for the interstellar 
absorption JVh. We can examine how this trial model fits the set of spectra by x 2 value, and 
we can improve the fit by changing the parameter values kT i: n gti , and Nr. This procedure is 
done with the XSPEC data analysis package as a simultaneous spectral fitting of N spectra. 
For some clusters, the interstellar absorption was poorly constrained. In such case, we fixed the 
absorption column to the Galactic value. We show the sample of the fitting result in Figure 1 
(Right). 

3. Temperature and Gas Density Profiles 

We determined the temperature and gas (electron) density profiles of all sample clusters. 
Figure 2 shows the temperature, gas density, pressure profiles of A2597. The pressure profiles 
are simply derived from the temperature and gas density with the equation of the state of ideal 
gas: P = n e kT. We attempted to model the temperature and density profiles with analytic 
functions. Note that the total mass profile can be calculated without employing such models, 
as shown in the next section. However, we investigated for another way to obtain the total 
mass profile, to which end such models are employed. We fitted the temperature profile with 
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exponential + constant model given by 

T(r) =T + T 1 exp(-r/r T ). (4) 

The fitting results are summarized in Table 3. Three parameters in the exponential + constant 
model were determined by the \ 2 fitting, but their error estimation was not trivial when the 
fitting was unacceptable. These large reduced \ 2 values are likely to be due to local fluctuations 
in the intrinsic temperature profiles or to unknown systematic errors in our analysis procedure. 
Note that we integrate the above model functions within one radius bin to obtain each model 
point. Thus, a coarse sampling is not the cause of the large x 2 values. In order to estimate 
conservative errors for the parameters, we assigned a systematic error to each data point in the 
temperature profile. The systematic error of the temperature is assumed to be the constant 
fraction of the measured temperature for all the data points, where the fraction is determined 
so as to obtain the reduced \ 2 °f unity in each temperature profile. Note that the total error is 
calculated to be the square root of the quadratic sum of the statistical and systematic errors. 
For some clusters, the original (i.e., before assignment of the systematic error) reduced y 2 is 
small enough to accept the fit. We assigned the systematic error only for the clusters for which 
the original fitting was rejected by the x 2 test with a significance level less than 1 %. We also 
fitted the gas density profile with NFW gas density model given by 

n g {r) = n e0 «p[-B(l - ^ffi^ )] (5) 

\r / r s ) 

(Makino, Sasaki, & Suto 1998). The fitting results are summarized in Table 4. As was the case 
for the temperature profiles, we adopted the systematic error for all the clusters. 

4. Mass Profiles 

Under the assumptions of hydrostatic equilibrium and spherical symmetry, we can obtain 

the total gravitating mass profile as a function of radius using 

dP„ d<p 
— = -»n g m p - 

GM(< r) 

= -/jn g m p — . (6) 

We derived this mass profile by two different methods in order to check its consistency. 
The first method employs the temperature and density profile models obtained in Section 
3. Substituting Equation 4 and Equation 5 for Equation 6, we obtain the mass profile in an 
analytic form. The second method does not employ the temperature and density profile models. 
Instead, the mass profile is derived by approximating Equation 6 as simple differences: 

1 r 2 AP(r) 

M(<r) J-— I_f±±H2. ( 7 ) 

\in g \r)m v G Ar 

We calculated AP(r)/Ar as 
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AP(r) = P l+1 -P t 

Ar r i+ i - n 

where and r; are the pressure and radius of the zth shell. The radius r and gas density n g (r) 
are given by r = (r i+ i + rj)/2 and n g (r) = (n e j+i +n ei )/2, respectively 

Results of these two methods are compared in Figure 3 for A2597. The plots also show 
the lcr confidence levels for the analytic mass profile, which is derived by considering the errors 
of the parameters describing the temperature and density profile models. The mass profiles 
derived by the two different methods were consistent in most of the cases. In the case that the 
pressure of the outer shell is larger than that of the inner shell (Pi+i > Pi), the mass profile 
shows the negative value at that point : M < 0. This is likely caused by the local temperature 
fluctuation intrinsic at some radius of the clusters or by systematic errors in our analysis. We 
excluded these unphysical points in the analysis, though such points are only seen in NGC5044 
and Centaurs. Of the two derivation methods, the first one using the modeled temperature 
and density profiles is easy to handle, but it involves sacrificing one important point of the 
deprojection analysis; i.e., no reliance on any particular profile models. On the other hand, the 
second method is more straightforward, but it suffers large error owing to local fluctuations in 
the temperature and density profiles. 

4-1- Scaling of Mass Profiles 

In Figure 4, we show the analytic mass profiles of 23 sample clusters in one plot, illustrat- 
ing the scatters among them. CDM simulations predict that the density profiles of dark matter 
are universal in form across a wide range of mass scales (Navarro, Frenk, & White 1995; Navarro, 
Frenk, & White 1996). We scaled our analytic mass profiles with r2oo and M200, where r2oo is 
the radius within which the mean halo density is 200 times the critical density of the universe, 
and M200 is the total mass enclosed within r2oo- As shown by Navarro, Frenk, & White (1995), 
and Navarro, Frenk, & White (1996), clusters of different mass are expected to show similar 
structures when scaled to such a characteristic radius and mass. For the calculation of r 2 oo, we 
used the relation obtained from the numerical simulation by Evrard, Metzler, & Navarro (1996): 

r 200 = 3.690 (T/10keV) a5 (1 + z)' 1 - 5 [Mpc], (9) 
where T is the spatially averaged temperature, and z is the redshift. M200 is calculated by 

M 200 = ^7r(200p crit (^)) r 2 3 00 . (10) 

where p CV it{z) = 3H(z) 2 /8ttG is the critical density of the universe at a redshift z. We show the 
scaled mass profiles in Figure 5. On a large scale (r > O.l^oo), the scaled mass profiles agree 
with each other better than did the original mass profiles, except in the case of one deviant 
profile of A401. This findings suggest that the mass profiles have a similar form on a large scale; 
in other words, the scaling with r2oo and M200 is effective at least on this scale. The standard 
deviation of the mass profiles is 41 % at 200 kpc for the original mass profiles, and that for 
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the scaled profiles is 21 % at 0.1r2oo ; which corresponds to about 160-300 kpc. In contrast, the 
standard deviations on the small scale (r < 0.1r2oo) are not significantly different: 55 % at 20 
kpc for the original mass profiles, and 60 % at 0.01r 2 oo for the scaled mass profiles. 

When the density profile of dark matter is described with the power-law expression 
p(r) = po( r / r o) _c \ the mass integrated over the volume is described by 

M(< r) = f 4irp(r')r' 2 dr' = ^MA(L)3-a. (n) 



3 — a ?~o ' 

Therefore, the smaller the value of a is (a — > 0), the steeper the mass profile is. We overlaid 
the M oc r 1 ' 5 (a = 1.5), M ocr 2 (a = 1) and M oc r 3 (a = 0) lines on the scaled mass profiles in 
Figure 5. It was found that the slope a was in the range of to 1.5, and it was flatter (smaller) 
on the small scale. The slope a at the cluster center is quantitatively examined in Section 4.2 

4-2. Inner Slope of Dark Matter Distribution 

The shape of the dark matter distribution near the center of a cluster is sensitive to 
the theoretical models adopted. In this section, we focus on the observed shape of the total 
mass distribution in terms of the slope of the density profile at the inner part of a cluster. The 
inner slope of the density profile is obtained by fitting the total mass profile we obtained with 
a model mass profile calculated from an assumed density profile. Although we employed the 
King and NFW profiles as density profile models in the previous section, a more generalized 
form is used in this section. That form was 

^( r ) = ( r / rs )«(l + ( r / rs ))(3-a)' ( 12 ) 

where p is the central density, r s is its scale radius, and a is the asymptotic slope of the 
profile at small radii. This form of density profile requires numerical integration to derive the 
integrated mass profile. The asymptotic slope a in Equation 12 can be used as the inner slope. 
However, it was found that the asymptotic slope a and the scale radius r s are coupled strongly, 
and therefore difficult to determine independently. Thus, we focus on the slope of the density 
profile at a finite radius and used it as the inner slope. The slope at a radius r , ao( r o) is given 
by 

Using Equation 12, we get 

/ \ Q + (r /r s )) 
(l + (r /r s )) 

We employ a (^o) instead of a, in addition to p an d r s , as a free parameter of the fitting. We 
fixed r to 0.02r 2 o(b which corresponds to about 40 kpc. The choice of r is not trivial, but we 
fix this value so that the radius is appropriate for a comparison of the theoretical models and is 
covered by observed data points in the mass profile. Nevertheless, A401, A644, and A963, due 
to a lack of data points within 0.02r 2 ooj were discarded from the following analysis. Since the 
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mass profiles of some clusters are nearly power-law, we could not determine r s . In such case, 
we fixed r s to 1 Mpc, which is the radius that there is no data for all clusters. 

The results of the fitting are summarized in Table 5, and the best-fit values and errors 
of the inner slope ao( r o) are plotted in Figure 6. We also show the total mass profiles with the 
best-fit models for the general form of the density profile in Figure 7. The inner slope «o spans 
a wide range with 0<ao< 2.3. We found that the 90% upper bound of a was lower than unity 
for 6/20 (-41 %) clusters (A2052, A2597, A478, PKS0745-191, ZW3146, and 2A0335+096), 
suggesting that the dark matter distribution in a significant fraction of clusters was flatter than 
that in CDM halo models such as the NFW profile or the Moore profile. 

5. Examination of the Systematic Effects 

We have demonstrated that the inner slope «o shows a large scatter and is less than 
unity for 30% of the clusters in our sample. Before discussing these results in greater detail, 
we would like to address their validity and the systematic effects that may affect the measured 
inner slope from various points of view. 

5.1. Center Position 

We can define three types of positions as the center of a cluster: (1) the X-ray emission 
peak of cluster hot gas; (2) the X-ray centroid; and (3) the position of the Brightest Cluster 
Galaxy (BCG). We chose the X-ray emission peak as the center of the annuli to extract spectra. 
Since these three positions are slightly different, the selection of the center position may affect 
our results. 

We first examined the position of the BCG in our sample clusters. The BCG is defined 
as the brightest galaxies among the member galaxies of a cluster. Since the BCGs are usually 
located at the center of a cluster and have velocities very near the mean velocity of galaxies in 
the cluster, they are considered to sit at the bottom of the cluster gravitational potential well. 
The positions of BCGs were taken from the NASA/IPAC Extragalactic Database (NED) 3 . We 
calculated the offset between the X-ray peak and the position of the BCG. When the offset was 
larger than or comparable to the radius of the innermost annulus, the temperature or density 
profiles will be affected by the selection of the center. In Figure 8, the offsets of the BCG are 
plotted against the radii of the innermost annulus. The offsets of the BCG are smaller than 
the radius of the innermost annulus (~ 30 % at the maximum) except for 2A0335, suggesting 
that the difference in cluster center between these two definitions did not affect our results 
significantly. 

We next examined the offset between the X-ray peak and the X-ray centroid. To derive 
the X-ray centroid, we used the X-ray images in which point sources were removed. We replaced 
each embedded source with the local diffuse X-ray emission surrounding the source by using 

3 http://nedwww.ipac.caltech.edu 
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the CIAO task dmfilth. This task replaces the counts within the source regions with the values 
sampled from the background regions. As the ACIS CCDs do not cover the entire X-ray 
emission for most clusters. We thus define the X-ray centroid within a circular region that is 
centered on the X-ray peak and has a radius that is tangent to the detector edge. The offsets 
between the X-ray peak and the X-ray centroid are plotted against the radii of the innermost 
annulus in Figure 9. The offsets of the X-ray centroid are slightly larger than the offsets of the 
BCG but are smaller than the radii of the innermost annulus. The offset is at most 74 % (for 
A644) of the radius of the innermost annulus, suggesting that the difference in these definitions 
of a cluster center does not significantly affect the results. 

5.2. Spherical Symmetry 

When we selected our sample, we excluded some clusters which are not spherical sym- 
metric in appearance. To quantify the spherical symmetry, we measured the ellipticity (e) and 
the position angle (PA) from the projected X-ray image. We used an iterative moment tech- 
nique derived from the treatment of the dispersion ellipse of the bivariate normal frequency 
function of position vectors used by Carter & Metcalfe (1980). We first calculated the mo- 
ments of the observed X-ray images. From an image of P pixels having rij counts in pixel i, we 
computed the moment 
1 p 

Vmn = -T7^2 , n i (xi-x) m (y i -y) n (m,n<2), (15) 
^ v i=i 

where N = J2[=i n i, an d (x,y) is the centroid. Then ellipticity e is 

and the position angle of the major axis measured north through east in celestial coordinates 
is 

PA = tan- 1 (-/H— ) + * (17) 

1\ + — /i 2 * 

where A± (A + > A_) are the positive roots of the quadratic equation 

^ 20_A2 ** u , =0. (18) 
/in A*02 - ^ 

As in the case of the determination of the X-ray centroids, we used the X-ray images from 
which point sources were removed and replaced with local diffuse X-ray emission surrounding 
the source. We also employed the ROSAT images to determine the ellipticities and position 
angles in outer regions of clusters. From the ROSAT image, the point sources were removed 
but the holes were not replaced with the background. We consider that if the ellipticity of a 
cluster affects the inner slope ao, the ellipticity e and the inner slope will show some kind 
of relation. However, the plot in Figure 10 shows no correlation between the ellipticity e at 
r = 0.05r2oo and the inner slope «o (the correlation coefficient is -0.31). We therefore conclude 
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that there was no evidence that deviation from spherical symmetry affects the results on the 
inner slope «o, although we cannot obtain no evidence regarding symmetry along the line of 
sight. 

5.3. Central Structure 

Recent Chandra observations have revealed remarkable structures in the hot gas of the 
central region of some cooling flow clusters. These structures may have observably affected the 
temperature and gas density in our measurement. Furthermore, if these structures indicate 
a break of hydrostatic equilibrium, they may systematically have affected the mass profile we 
obtained. For 10 clusters in our sample, the presence of central X-ray structures such as a 
cavity, hole, or plume has been reported with Chandra observations in the literature (Table 6). 

In Figure 11, we show the inner slope «o again, indicating 10 clusters for which the 
central structure has been found by open circles. It is found that the three clusters in which 
«o is as steep as 2 have central structures, and the range of cto becomes narrower if we neglect 
them. However, the distribution of ao from to 1.2 is similar for clusters with and those 
without central structures, though this is difficult to conclude quantitatively. Note that some 
of the central structures, such as those in A133, A2597, and MKW3S are small enough to be 
removed from the analysis. In order to evaluate the observational effect of these structures, we 
removed the region of these structures from those clusters and confirmed that the mass profiles 
were not significantly affected. 

6. Relations between Inner Slope and Other Observational Parameters 

As shown in Section 4.2, the inner slope «o of the density profiles spans a wider range 
than that estimated from their errors. Even if we neglect the three clusters in which central 
structures might affect the results, «o ranges from to 1.2 and spreads toward a flatter side than 
expected based on CDM simulations. If this spread of the distribution is intrinsic, what is it 
that determines the inner slope of the density profile? In this section, we explore observational 
parameters that primarily determine the inner slope ao, by examining their correlations. 

6.1. Redshift vs. «o 

We first show the relation between the redshift and the inner slope «o in Figure 12. The 
correlation coefficient is -0.25 for this relation. Note that the range of redshifts in our sample 
may not be sufficient to investigate the effect of cosmological evolution. 

6.2. Temperature vs. ciq 

We plot the spatially averaged temperature fcT average against the inner slope a in Figure 
13. The correlation coefficient is 0.01. Simple arguments based on virial theorem suggest that 
the mass of a cluster is simply related to the cluster temperature as Mix T 3//2 . This relation 
implies that the inner slope ao is not related to the scale of the cluster. 
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6.3. Gas Fraction vs. «o 

We next investigated the relation between the gas fraction and the inner slope eta. The 
gas fraction is the ratio of the hot gas mass to the total mass, and is defined as a function of 
radius. The integrated gas mass profile M gas (< r) is given by 



where n gas (r) is the total number density of electrons and ions, fi(= 0.6) is the mean molecular 
weight, and m p is the proton mass. The gas fraction is defined as 



In Figure 14, we present the profiles of total mass, gas mass, and gas fraction for 23 
sample clusters. It was found that the gas fraction increases toward the center for some clusters. 
In Figure 15, the gas fractions at the radius of r = 0.05r 20 o are plotted against the inner slope a^. 
A negative correlation was observed with the correlation coefficient of —0.51, for which case, at 
a significance level of about 3%, the hypothesis of no correlation is rejected. This correlation 
might be a kind of artifact in the analysis, since we derive both the total mass profile and the 
gas mass profile from the same gas density profile and gas temperature profile. However, it is 
unlikely that the observed correlation is due to correlated errors between the two parameters, 
considering the size of the errors. 

In order to examine this correlation is artifact or not, we took M 2 oo, which is deter- 
mined solely from the gas temperature and redshift, instead of the integrated mass profile. We 
redefined the gas fraction as 



We plot /g as at the radius of r = 0.05r 2 oo against the inner slope a in Figure 16. Although 
the correlation coefficient of -0.40 is smaller than that for / gas , the no correlation hypothesis is 
rejected at a significance level of less than 10%. 

Correlations between the inner slope ao and the gas fraction according indicate that 
gas-rich clusters in the central region tend to have a flat core: a < 1. 

7. Summary and Discussion 

We have analyzed the Chandra data of 23 clusters of galaxies in order to investigate 
central mass distribution. The high spatial resolution imaging spectroscopy of Chandra and 
a new deprojection technique enable us to measure the temperatures and gas densities in the 
very central region of sample clusters without assuming any particular models. Under the 
assumptions of hydrostatic equilibrium and spherical symmetry, we obtained the deprojected 




(19) 




(20) 




(21) 
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mass profiles. Our major results are as follows. 

1. The mass profiles scaled with r2oo and M200 agree each other on the large scale r > 0.1r2oo- 
In contrast, the central (r < 0.1r 2 oo) mass profiles show a large scatter. 

2. The inner slope of the density profile was derived by fitting the mass profile with a 
general form of dark matter density profile for 20 clusters. The values of ao span a wide 
range (from to 2.3). For 6 out of 20 clusters, «o are lower than unity at a 90 % confidence 
level. 

3. We investigated several features that might influence the results of the inner slope, in- 
cluding center position, ellipticity, and central structure of a cluster. We found that the 
systematic effects of these features are not significant except in the case of the central 
structure, which may broaden the distribution of the inner slope ao- However, even if we 
excluded the clusters in which central structures were seen, the inner slope «o distributes 
widely in the range from to 1.2. 

4. We examined the relationships between the inner slope ao and other observational param- 
eters. Although redshift, averaged temperature, and variation in the temperature profile 
are not correlated with the inner slope «o, gas fraction near the center of a cluster has a 
negative correlation with ao- 

CDM simulations predict that the inner slope a is in the range 1 < a < 2. Therefore, 
our results are inconsistent with the CDM simulations. Our observations provide flatter slopes, 
at least for some clusters, than those expected from the CDM simulations. This is true even if 
we neglect the clusters showing central structures (0 < a < 1.2). 

A similar claim has been presented by authors who are investigating the rotation curve 
of galaxies and clusters. Firmani et al. (2001) examined the observed rotation curves of dwarf 
and low surface brightness (LSB) galaxies, and two clusters of galaxies, and found that all 
of those objects have soft cores: a < 1. Swaters, Madore, van den Bosch, & Balcells (2003) 
observed the rotation curve of 15 dwarf and LSB galaxies and found inner slopes in the range 
of ^ a < 1 in the majority. This inconsistency between observations and simulations in terms 
of the dark matter distribution at the center of galaxies or clusters is called the core problem. 
Our results indicate that the core problem exists in a significant fraction of clusters observed 
through X-ray observations. 

We also found a negative correlation between the inner slope and the gas fraction. In 
the central region of clusters, the baryonic components such as hot gas and stars in a galaxy are 
not negligible. These baryonic components are usually considered to follow the gravitational 
potential, which is predominantly determined by dark matter distribution. However, our results 
imply that these baryonic components affect the dark matter distribution in the central region 
of some clusters. 

El-Zant, Shlosman, & Hoffman (2001) argued that the core problem can be resolved 
within the framework of the standard CDM model by considering the dynamical friction (DF) 
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between dark matter and gas. El-zant et al. assumed that the gas is not initially smoothly 
distributed but is, rather, concentrated in clumps. Such gas clumps move through smooth 
dark matter particles, lose energy to the central dark halo, and heat it up. This leads to the 
puffing up of the central regions and to the flattening of the density profile. Monte Carlo 
simulation by El-zant et al. successfully reproduced the observed flat density profile. In the 
DF model, the baryonic (gas) component is more centrally concentrated than dark matter, 
because gas gives its energy to dark matter and shrinks toward the center. So this model 
accounts for our observational result that gas rich clusters in the central region tend to have 
a flat core. Note that another numerical simulations including baryonic components show 
the cooling baryons adiabatically compress the dark matter and make a steeper central density 
profile (e.g., Blumenthal, Faber, Flores, & Primack 1986). However, this result is not consistent 
with our results because the gas rich cluster in the central region is expected to have steeper 
profile if the cooling baryons adiabatically compress the dark matter. 

Although the DF model is one of many interpretations proposed to account for the core 
problem, our results indicate that the numerical simulation need to reconsider the treatment 
of the baryonic components in galaxies or galaxy clusters. 

We thank all the staff members involved in the Chandra project. We acknowledge Jone 
Arabadjis and Mark Bautz for their support on the deprojection technique and their useful 
comments. 
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Table 1. Observation log of the sample clusters. 



ID 


Cluster 


Obs. 


Ra (deg)° 


Dec (deg)° 


Observation date 


Exp. 


Screened 6 


ACIS C 






ID 








(ks) 


Exp. (ks) 


chip 


1 


A1060 


2220 


159.073 


-27.569 


2001-06-04 04:43:23 


31.9 


30.0 


I 


2 


A133 


2203 


15.689 


-21.882 


2000-10-13 22:27:02 


35.5 


34.5 


S 


3 


A1795 


493 


207.205 


26.608 


2000-03-21 07:54:49 


19.6 


19.6 


S 






3666 


207.204 


26.576 


2002-06-10 16:21:19 


14.4 


13.8 


s 






494 


207.236 


26.607 


1999-12-20 05:00:57 


19.5 


17.6 


s 


4 


A1835 


495 


210.272 


2.895 


1999-12-11 16:48:33 


19.5 


19.5 


s 






496 


210.222 


2.867 


2000-04-29 06:55:44 


10.7 


10.7 


s 


5 


A2029 


891 


227.725 


5.764 


2000-04-12 06:38:56 


19.8 


19.8 


s 


6 


A2052 


890 


229.182 


7.012 


2000-09-03 06:01:22 


36.8 


36.8 


s 


7 


A2199 


498 


247.188 


39.553 


1999-12-11 10:47:37 


18.9 


18.9 


s 






497 


247.135 


39.560 


2000-05-13 17:36:15 


19.5 


17.9 


s 


8 


A2204 


499 


248.185 


5.557 


2000-07-29 02:49:42 


10.1 


10.1 


s 


9 


A2597 


922 


351.337 


-12.135 


2000-07-28 05:13:47 


39.4 


25.1 


s 


10 


A401 


518 


44.727 


13.579 


1999-09-17 21:35:26 


18.0 


18.0 


I 






2309 


44.732 


13.461 


2000-11-03 19:10:36 


11.6 


11.6 


I 


11 


A478 


1669 


63.362 


10.436 


2001-01-27 03:28:03 


42.4 


42.4 


I 


12 


A644 


2211 


124.329 


-7.543 


2001-03-26 00:27:49 


29.7 


29.2 


s 


13 


A85 


904 


10.442 


-9.374 


2000-08-19 07:06:52 


38.4 


38.4 


I 


14 


A963 


903 


154.284 


39.063 


2000-10-11 00:01:18 


36.3 


36.3 


s 


15 


AWM7 


908 


43.665 


41.664 


2000-08-19 18:30:01 


47.9 


47.9 


I 


16 


Centaurus 


504 


192.207 


-41.334 


2000-05-22 00:33:17 


31.7 


28.8 


s 






505 


192.199 


-41.334 


2000-06-08 00:51:50 


10.0 


10.0 


s 


17 


Hydra A 


575 


139.527 


-12.091 


1999-10-30 07:29:02 


23.8 


23.8 


I 






576 


139.527 


-12.091 


1999-11-02 11:31:54 


19.5 


19.5 


s 


18 


MKW3S 


900 


230.488 


7.757 


2000-04-03 12:26:13 


57.3 


57.3 


I 


19 


NGC5044 


798 


198.859 


-16.378 


2000-03-19 15:42:42 


20.5 


19.8 


s 


20 


PKS0745-191 


2427 


116.860 


-19.306 


2001-06-16 05:32:52 


17.9 


17.9 


s 






508 


116.870 


-19.277 


2000-08-28 22:15:31 


28.0 


4.6 


s 


21 


Sersicl59-03 


1668 


348.515 


-42.713 


2001-08-13 16:41:20 


9.9 


9.9 


s 


22 


ZW3146 


909 


155.905 


4.166 


2000-05-10 03:20:25 


46.0 


46.0 


I 


23 


2A0335+096 


919 


54.666 


10.008 


2000-09-06 00:03:13 


19.7 


14.1 


s 



a Nominal pointing position of the observation in Equinox 2000.0 
b Exposure time after lightcurve screening (see §4.2) 
c Detector on the aim point 



Table 2. Properties of the sample clusters. We show the redshifts, hydrogen column densities of the galactic absorption, 
temperatures, and X-ray fluxes of the sample clusters. The temperatures are referred to Reiprich & Bohringcr (2002), Ota 
(2000), and Allen et al. (1996). 

Cluster redshift N£ kT f x b Ref. c 







[10 20 cm- 2 ] 


[keV] 


[10 11 ergs s 1 cm 2 ] 




A1060 


0.0126 


4.79 


o O/I+0 06 
°-^-0.06 


9.95 


R 


A133 
A1795 


0.0570 
0.0622 


1.55 
1.20 


3.80t»;°° 
7.80!!;°° 


2.12 
6.27 


R 
R 


A1835 


0.2530 


2.30 


7.42™ 


1.47 





A2029 
A2052 


0.0780 
0.0345 


3.07 
2.78 


q iQ+i.oo 

q nQ+ 04 
°- uo -0.04 


6.94 
4.71 


R 
R 


A2199 
A2204 


0.0310 
0.1511 


0.87 
5.66 


4 i n+0-08 

7 Qi +0.25 
' - zl -0.25 


10.64 
2.75 


R 
R 


A2597 


0.0822 


2.50 


4.4018$ 


2.21 


R 


A401 


0.0748 


10.3 


8-00l°;l° 


5.28 


R 


A478 


0.0881 


14.8 


8-40l°io 


5.15 


R 


A644 


0.0704 


6.95 


7 qo+o- 80 

' " yu -0.80 


4.02 


R 


A85 


0.0557 


3.37 


6.90l°;f 


7.43 


R 


A963 


0.2057 


1.40 


6.831°;^ 


0.59 





AWM7 


0.0172 


9.91 


q ? r+0.09 

°- ' d -0.09 


1.58 


R 


Ccntraurus 


0.0110 


8.07 


3.681°;°^ 


27.19 


R 


Hydra A 


0.0538 


4.90 


4 30+°- 4 ° 


4.78 


R 


MKW3S 


0.0450 


3.04 


q 7n +0.20 
°- ,u -0.20 


3.30 


R 


NGC5044 


0.0089 


5.03 


1 07+°- 01 
i - u '-o.oi 


5.51 


R 


PKS0745-191 


0.1028 


40.7 


7 qi+0.11 

' - zl -o.n 


2.44 


R 


Scrsicl59-03 


0.0580 


1.76 


3 00+ 1 20 


2.49 


R 


ZW3146 


0.2906 


2.94 


6 10+ 030 

o - lu -0 .30 


0.66 


A 


2A0335+096 


0.0349 


17.6 


3 01 +°- 07 


9.16 


R 



a Hydrogen column density of the galactic absorption. 

b X-ray flux in units of 10 — 11 ergs s _1 cm -2 . The energy bands are 0.1-2.4 keV for R and 2-10 keV for O and A. 
c References R:Reiprich & Bohringcr (2002), 0:Ota (2000), and A:Allen et al. (1996). 



16 



Table 3. Fitting results of the temperature profiles with the exponential + constant model. 



Cluster 


T [keV] a 


Ti [keV] a 


rf. [kpc] 


xLt/dof b 


Systematic 
error c 


A1060 


3.16 ± 0.22 


-0.00 ± 0.80 


850.0 ± 918.5 


126.4/4 


0.183 


A133 


4.52 ± 0.36 


-3.45 ± 0.33 


90.9 ± 28.5 


20.7/4 


0.082 


A1795 


6.03 ± 0.51 


-3.36 ± 0.57 


151.4 ± 63.1 


105.1/5 


0.099 


A1835 


10.10 ± 0.69 


-7.60 ± 0.56 


103.4 ± 28.2 


10.2/4 




A2029 


9.89 ± 0.30 


-5.73 ± 0.34 


144.4 ± 22.7 


3.9/4 




A2052 


3.20 ± 0.27 


-5.67 ± 8.39 


14.9 ± 22.4 


24.8/2 


0.734 


A2199 


4.92 ± 0.32 


-3.34 ± 0.41 


57.2 ± 22.2 


97.3/7 


0.093 


A2204 


8.45 ± 1.29 


-7.21 ± 1.80 


58.1 ± 34.0 


26.1/3 


0.235 


A2597 


6.02 ± 0.39 


-4.33 ± 0.35 


162.7 ± 29.1 


6.6/4 




A401 


7.71 ± 0.62 


0.00 ± 29.09 


12.9 ± 0.1 


19.0/2 


0.170 


A478 


7.03 ± 0.35 


-5.52 ± 0.78 


48.4 ± 20.0 


33.6/6 


0.076 


A644 


6.47 ± 0.39 


-28.72 ± 28.72 


14.9 ± 27.3 


27.6/3 


0.128 


A85 


6.44 ± 0.27 


-4.56 ± 0.46 


78.4 ± 21.6 


26.0/5 


0.065 


A963 


6.22 ± 0.25 


-29.78 ± 23.80 


13.3 ± 6.1 


1.5/2 




AWM7 


3.78 ± 0.10 


-29.60 ± 27.94 


7.2 ± 3.8 


52.8/4 


0.064 


Centaurus 


4.63 ± 0.90 


-3.92 ± 0.86 


78.0 ± 35.2 


1167.4/7 


0.215 


Hydra A 


3.80 ± 0.24 


-0.95 ± 0.28 


113.3 ± 95.0 


33.4/4 


0.064 


MKW3S 


3.71 ± 0.14 


-1.10 ± 0.59 


30.9 ± 23.2 


42.3/7 


0.078 


NGC5044 


1.65 ± 0.39 


-1.10 ± 0.37 


102.9 ± 54.5 


168.6/7 


0.065 


PKS0745-191 


11.16 ± 0.65 


-8.89 ± 0.57 


154.6 ± 25.7 


9.9/3 




Sersicl59-03 


2.69 ± 0.09 


-1.31 ± 0.13 


58.0 ± 15.7 


1.8/2 




ZW3146 


8.04 ± 0.31 


-6.86 ± 0.58 


57.2 ± 13.0 


0.6/2 




2A0335+096 


3.96 ± 0.40 


-3.18 ± 0.35 


97.1 ± 30.2 


60.6/4 


0.089 



a Errors are estimated by including a systematic error (see text). 
b Original reduced \ 2 without a systematic error (Xstat)- 

c Systematic error adopted so as to get the reduced x 2 value of unity. This error includes the fluctuations in the intrinsic 
temperature profiles or unknown systematic errors in our analysis procedure. 
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Table 4. Fitting results of the gas density profiles with the NFW gas density model. 



Cluster 


n c0 [1(T 2 cm- 3 ]" 


r s [kpc] a 




xLt/dof b 


Systematic 
error c 


A1060 


0.83±0.09 


145.1±48.2 


5.13±1.08 


565.3/4 





,078 


A133 


4.21±0.80 


91.6±21.9 


6.40±0.37 


1584.3/4 





,121 


A1795 


4.76±0.45 


252.8±37.6 


8.54±0.50 


6460.0/5 





,084 


A1835 


22.10±1.52 


112.1±9.3 


7.04±0.16 


257.3/4 


0, 


,051 


A2029 


5.84±0.08 


193.7±4.5 


7.11±0.06 


31.9/5 





Oil 


A2052 


4.74±0.24 


64.5±5.3 


5.72±0.13 


144.1/2 





,029 


A2199 


3.00±0.26 


235.1±40.6 


8.52±0.66 


5112.0/7 





,079 


A2204 


21.60±3.10 


89.3±14.4 


7.15±0.28 


557.3/3 





,10 


A2597 


8.11±0.62 


157.4±17.6 


8.94±0.38 


537.9/4 





,071 


A401 


0.71±0.03 


2842.0±664.9 


19.05±10.64 


97.6/2 





,035 


A478 


6.46±1.25 


338.3±118.0 


9.99±1.53 


13901.8/6 





,155 


A644 


1.63±0.17 


493.6±125.8 


8.30±1.05 


972.4/3 





,085 


A85 


3.86±0.48 


122.5±21.6 


5.76±0.26 


2246.4/5 





,081 


A963 


2.76±0.13 


477.2±58.5 


9.07±0.57 


39.3/2 





,033 


AWM7 


1.26±0.14 


140.7±17.7 


5.00±0.26 


2899.1/4 





,083 


Centaurus 


7.79±1.75 


24.7±5.5 


5.72±0.23 


6499.1/7 





,143 


Hydra A 


7.03±1.53 


87.9±22.7 


6.67±0.43 


6555.0/4 





,151 


MKW3S 


3.23±0.33 


106.7±20.5 


5.95±0.41 


608.8/7 





,083 


NGC5044 


3.98±0.95 


34.8±15.5 


6.22±0.95 


3186.8/7 





163 


PKS0745-191 


12.00±0.38 


145.6±8.1 


7.39±0.16 


72.3/3 





,022 


Scrsicl59-03 


5.14±0.40 


240.5±54.8 


11.08±1.46 


101.5/2 





,057 


ZW3146 


14.30dzl.99 


258.2±63.0 


9.74±1.00 


473.6/2 





096 


2A0335+096 


9.23±1.30 


90.5±17.3 


7.41±0.47 


1887.5/4 





111 



a Errors are estimated by including a systematic error (see text). 
b Reduced x 2 without a systematic error (Xstat)- 

c Systematic error adopted so as to make the reduced \ 2 value unity. 
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Table 5. Fitting results of the mass profile with the general form of density profile given by Equation 12. The errors are 90 
% confidence intervals. 



Cluster r s [kpc] a a (0.02r 2 oo) a xLt/dof 6 Systematic 



error c 


A1060 


6.0±4.2 


1.17±0.38 


120.3/3 


0.224 


A133 


1000 (fixed) 


1.32±0.38 


31.8/3 


0.398 


A1795 


1000 (fixed) 


1.29±0.43 


194.2/4 


0.482 


A1835 


51.1±200.1 


0.82±0.47 


13.9/3 


0.278 


A2029 


760.0±1733.0 


1.15±0.47 


7.0/3 




A2052 


3.4±8.1 


-0.14±0.96 


176.4/1 


0.524 


A2199 


163.3±8.2 


0.64±0.46 


36.8/6 


0.440 


A2204 


1000 (fixed) 


1.13±0.56 


19.0/2 


0.525 


A2597 


21.8±76.6 


0.52±0.33 


20.9/3 


0.192 


A478 


3.1±3.5 


0.17±0.42 


23.8/4 


0.308 


A85 


1000 (fixed) 


0.94±0.37 


95.6/4 


0.365 


AWM7 


7.1±8.1 


0.48±0.39 


70.7/3 


0.290 


Centaums 


10.0±20.8 


2.28±0.46 


549.7/4 


0.640 


Hydra A 


1000 (fixed) 


1.85±0.55 


237.3/3 


0.556 


MKW3S 


86.4±101.8 


1.22±0.18 


17.1/4 




NGC5044 


1000 (fixed) 


1.76±0.21 


32.8/5 


0.217 


PKS0745-191 


17.4±6.6 


0.68±0.18 


2.3/2 




Sersicl59-03 


64.5±25.0 


0.68±0.67 


15.1/1 


0.378 


ZW3146 


69.4±9.0 


0.19±0.20 


0.3/1 




2A0335+096 


15.4±35.3 


0.56±0.10 


11.9/3 





a Errors are estimated by including a systematic error (see text). 

b Reduced \ 2 before including a systematic error (xftat)' 

c Systematic error adopted so as to make the reduced x 2 value unity. 
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Table 6. Remarkable structures in the central region of the sample clusters in literatures. 



Cluster 


Structure 


Reference 


A133 


tongue 


Fujita et al. (2002) 


A1795 


filament 


Fabian ct al. (2001) 


A2052 


holes 


Blanton, Sarazin, McNamara, & Wise (2001) 


A2199 


depression 


Johnstone, Allen, Fabian, & Sanders (2002) 


A2597 


cavities 


McNamara et al. (2001) 


Centaurus 


plume 


Sanders & Fabian (2002) 


Hydra A 


depression 


McNamara et al. (2000) 


NGC5044 


hole 


Buote, Lewis, Brighenti, & Mathews (2003) 


MKW3S 


filament & depression 


Mazzotta et al. (2002) 


2A0335+096 


cavities & blobs 


Mazzotta, Edge, & Markevitch (2003) 
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A1835 # 




Fig. 1. (Left) Schematic view of the deprojection analysis. The left spheres represent X-ray emissions 
from N spherical shells, and the right circles represent observed projected luminosities which are the 
integrations of emissivities along the line of sight. (Right) An example of the spectrum fitting. In this 
example (A1835), we fitted the N = 7 spectra for two data sets simultaneously. Observation IDs of each 
data set are 495 and 496. 



A2597 

o p 1 1 1 1 — i — i — i — ^ 




Fig. 2. The temperature, electron density, and pressure profiles of A2597. The solid lines represent best-fit 
models with analytical functions (Equation 4 and Equation 5). The pressure is simply derived from the 
temperature and gas density with the equation of state of the ideal gas: P = n e kT. 
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Fig. 3. The mass profiles of A2597. The solid curves represent the mass profile derived from the best-fit 
parameters of Equation 4 and Equation 5. The dashed lines are a confidence level of 68 %. The discrete 
data with error bars represent the mass profile calculated using approximate expression given by Equation 
7. 




10 100 
Radius [kpc] 



Fig. 4. Mass profiles of 23 sample clusters obtained from the temperature and density profile models. 



S 9 
a' ° 

3 



o 




0-01 0.1 

r ' /r 200 

Fig. 5. Mass profiles scaled by A/200 and r2oo- For the calculation of r2oo and M200, we used the relation 
obtained from the numerical simulation by Evrard, Metzler, & Navarro (1996). The dashed lines represent 
Mocr 1 - 5 (a = 1.5), Mocr 2 (a = 1) and Mocr 3 (a = 0), respectively. 
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5 10 15 20 

§ of cluster 

Fig. 6. Values of the Inner slope ao at the radius of 0.02r2oo for 20 clusters in Table 5. A401, A644, and 
A963 were removed due to a lack of data points within 0.02r2oo- The horizontal dashed lines represent 
a = 1.5 (Moore), a = 1.0 (NFW), and a = 0.0 (King). Error bars are shown at 90 % confidence level. 
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Fig. 7. The best-fit mass models using the general form of the density profile. The horizontal dashed 
lines represent the radius of ro = 0.02r2oo- The dash-dotted lines represent the best-fit slope (M oc r 3 ~ a ° 
at the radius r . A401, A644, and A963 are removed because of the lack of the data within 0.02r 2 oo- 
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Fig. 8. Radius of the innermost annulus vs. offset between the X-ray peak and the BCG. The dashed 
line represents 50 % of the radius of the innermost annulus. 




Fig. 9. Radius of the innermost annulus vs. offset between the X-ray peak and the X-ray centroid. The 
dashed line represents 50 % of the radius of the innermost annulus. 
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Fig. 10. Ellipticity e vs. inner slope ao. The errors are la (68.3 %) confidence level. The correlation 
coefficient is —0.31 for this relation. 
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Fig. 11. Inner slope ao of mass profiles. The circles indicate 10 clusters in the literatures to have the 
central structure (see Table 6). 



0.15 
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Fig. 12. Redshift vs. inner slope cto. The correlation coefficient is -0.25 for this relation. 




Fig. 13. Temperature vs. inner slope «o- The correlation coefficient is 0.01 for this relation. 
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Fig. 14. Total integrated mass, gas mass, and gas fraction profiles of the 23 sample clusters. The radius 
is scaled with t^qq- We show the total integrated mass and gas mass profiles in the upper panel, and the 
gas fraction profile in lower panel. 
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Fig. 15. Gas fraction at the radius of r = 0.05r2oo vs. inner slope «o- The correlation coefficient is — 0.5i 
for this relation. 
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Fig. 16. Gas fraction redefined in Equation 21 (/LJ vs. inner slope «o- The correlation coefficient is 
—0.40 for this relation. 
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